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Abstract 

Stationary and isotropic iteration stable random tessellations are considered, 
which can be constructed by a random process of cell division. The collection of 
maximal polytopes at a fixed time t within a convex window C M'^ is regarded 
and formulas for mean values, variances, as well as a characterization of certain co- 
variance measures are proved. The focus is on the case d >3, which is different from 
the planar one, treated separately in [13]. Moreover, a multivariate limit theorem 
for the vector of suitably rescaled intrinsic volumes is established, leading in each 
component - in sharp contrast to the situation in the plane ~ to a non-Gaussian 
limit. 
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1 Introduction 



Random tessellations stable with respect to iterations (STIT) have recently been intro- 
duced in [6] and [8] as a new model for random tessellation in and have quickly attracted 
considerable interest in modern stochastic geometry as well fitting the growing demand for 
non-trivial and flexible but mathematically tractable tessellation models. The STIT tessel- 
lations may be interpreted as outcome of a random cell division process, which makes them 
very attractive for applications, see for example [10]. Other potential applications include 
mathematical modelling of systems of cracks, joints or fissures in rock or the so-called 
craquelee of thin layers. A general approach to random cell division processes has recently 
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appeared in [2] and the construction of STIT tessellations can roughly be described as fol- 
lows. At first, we fix a compact and convex window W gM.'^ in which the construction is 
carried out. For simplicity we assume to be a polytope, as this implies that the resulting 
tessellation has polytopal cells with probability one. Next, an exponentially random life 
time is assigned to W, whereby the parameter of the distribution is given as a constant 
multiple of the integral-geometric mean width of W, see Section 2 for details. Upon expiry 
of this life time a W hitting uniform random hyperplane is chosen, is introduced in W and 
divides the window into two polytopal sub-cells W'^ and W~. The construction continues 
now recursively and independently in both of these sub-cells, where the newly introduced 
hyperplanes are always chopped-off by the boundaries of their mother-cells. The whole 
construction is continued until some deterministic time threshold t is reached. Regarded 
in time, the construction can by interpreted as a pure-jump Markov process on the space 
of tessellations of the window W. 

Our construction shares some common features with random fragmentation processes or 
branching Markov chains in the sense of [1]. The cells of the resulting tessellation within W 
can be regarded as particles in a suitable Pohsh space and the dynamics of the particles is 
non-interacting in a sense that different particles (cells of the tessellation) have independent 
evolutions, which is indeed the case in our construction. Moreover, whenever a particle 
dies it is replaced by exactly two new particles, namely the two newly generated sub-cells. 
In addition, the hfe times of the particles are exponentially distributed, as assumed in [1], 
but in general not independent in contrast to the fragmentation theory. 
In the recent paper [12] the authors have introduced a new technique for studying the 
geometric properties of STIT tessellations based on martingales and the general theory of 
martingale problems for pure-jump type Markov processes. In particular, with these new 
developments, the variance of the total surface area of a stationary and isotropic iteration 
stable random tessellation Y{t, W) in 1^ C M*^ (a random STIT tessellation) has been 
determined by integral-geometric means and the corresponding central limit theory has 
been established. Strikingly, as already signalled by results in the special two dimensional 
case [13], and as confirmed by the results of the present paper, it turns out that in a cer- 
tain rather strong sense the asymptotic behaviour of the surface area process dominates 
and fully determines the asymptotic geometry of the STIT tessellation. In particular the 
surface area variance is the basic second-order parameter of the tessellation and all second- 
order characteristics of functionals considered in our work can be reduced to it, likewise 
non-trivial functional limits in law exist for the STIT surface area process whereas the 
limits in law for other natural related processes arise as deterministic functionals of the 
corresponding surface process. Another crucial phenomenon arising in this context, as first 
noted in [12], is that the asymptotic theories for dimension d = 2 and d > 2 differ strongly 
in many important aspects. We have studied the planar case o? = 2 in the recent separate 
paper [13] and thus only consider d > 2 in the present one. 

The purpose of this paper is to establish a second-order and limit theory for integral- 
geometric characteristics of stationary and isotropic STIT tessellations in dimensions higher 
than two. More precisely, the characteristics stiidied in this work are the cumulative in- 
trinsic volumes of all orders for the collections of so-called maximal polytopes of Y{t, W). 
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These are in codimension 1 the cell-separating facets introduced during the random cell di- 
vision process described above, constituting the basic building blocks of a STIT tessellation 
as discussed in detail below. We shall provide explicit as well as asymptotic variance ex- 
pressions for these parameters of the random tessellation Y{t, Wr) as R tends to infinity for 
a sequence Wr — RW of expanding convex windows. Further, we will find the covariance 
measures for random lower-dimensional face measures generated by Y{t,W). Finally, we 
shall also give the corresponding convergence in law statements, obtaining non-Gaussian 
limits for the studied case c? > 2 as opposed to the classical Gaussian limits arising for 
d= 2, see [13]. 

The paper is organized as follows: In Section 2 we define STIT tessellations, specializ- 
ing to the stationary and isotropic set-up in the focus of this paper, and we discuss their 
basic properties that are needed in our arguments. Next, in Section 3 we calculate the vari- 
ances of the cumulative intrinsic volumes of all orders for STIT tessellations. This includes 
both, exact formulae and asymptotic analysis upon letting the window size grow to infinity. 
Further, in Section 4 we extend the second-order analysis to the level of lower-dimensional 
face measures induced by STIT tessellations, thus taking into account not only the numeric 
characteristics but also the spatial profile of the STIT face processes. Finally, in Section 
5 we develop the corresponding functional limit theory with non-Gaussian limit processes 
for d > 2. In order to keep the paper self-contained, we will recall important faces from 
[12] and sketch some of their proofs for the readers convenience. 

2 STIT Tessellations 

The purpose of this Section is to provide a short self-contained discussion of STIT tes- 
sellations in M.'^ as studied in this paper and to summarize their basic properties for easy 
reference. We will restrict to the stationary and isotropic case in Subsection 2.1 below, 
speciahzing to the scope of the paper. The general reference for these - by now classical - 
properties throughout this section is [8]. Next, in Subsection 2.2 we will discuss the mar- 
tingale tools developed in [12] and underlying our present theory. Finally, in Subsection 
2.3 we provide certain useful mean value relationships for intrinsic volumes in context of 
STIT tessellations. 

2.1 The MNW-Construction and Basic Properties 

We start with a compact and convex polytope W gW^ in which our construction is carried 
out and denote by A the standard isometry-invariant measure on the space T-L of (affine) 
hyperplanes in M.^ normalized so as to induce unit surface intensity on H (note that in 
our earlier papers [12] and [13] this measure has been denoted by Aj^o)- We call A the 
driving measure of the construction. Assign now to W an exponentially random lifetime 
with parameter A([Vr]), where 

[W]:^{H en-. Hnw 
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is the set of hyperplanes hitting W. Upon expiry of this random hfe time, a random hyper- 
plane is chosen according to the distribution A([1^])~^A(- fl [W]), is introduced in W and 
is chopped off by its boundary. This is, the window W sphts into the two polyhedral sub- 
cells and W~ that are separated by the introduced hyperplane piece. The construction 
continues now recursively and independently in and W~ and is stopped if some pre- 
viously fixed deterministic time threshold i > is reached. Our assumptions ensure that 
the cells of the tessellation constructed until time t > are convex polyhedra in W with 
probability one. They are denoted by Cells(y(t, W^)) and we denote by Y{t,W) the ran- 
dom closed set in W that consists of the union of cell-boundaries of cells constructed until 
time t, see Figure 1. The construction of Y{t, W) is referred to as the MNW- construction 
after the names of its inventors - Mecke, Nagel and Weiss - and the tessellation Y{t, W) 
itself is called a random STIT tessellation. The abbreviation STIT comes from the cru- 
cial property enjoyed by Y{t, W)^ namely that of being stable under iteration, for which 
we refer to [6] or [8]. We call the cell separating (d — l)-dimensional faces the {d — 1)- 
dimensional maximal polytopes and denote the collection of all such polytopes of Y{t, W) 
by MaxPolytopeS(^_i(y(i, 1^)). Moreover, we introduce the set MaxPolytopeSj(y(i, 1^)) 
of j-dimensional faces of {d — l)-dimensional maximal polytopes of Y{t, W) by 

MaxPolytopes^(F(t,H/)) = |J FaceSj(/), 

/GMaxPolytopeSrf_i(y(t,VK)) 

where by FaceSj(/) we mean the set of all j-dimensional faces of the (d — 1) -dimensional 
polytope /, < j < d — 1. The maximal polytopes are often referred to as I-polytopes 
- this terminology originates from the historically first considered particular case d — 2 
where the maximal polytopes (which are just line segments in this case) assumed shapes 
similar to the literal I. 

It is an important observation that the spatio-temporal construction of Y{t, W) satisfies 
the Markov property in time parameter t, which means that 

Y(t + s,W)^ Y(t, W) ffl Y(s, W), 

where ffl denotes the operation of iteration of tessellations, see [8], [6] or [12]. To make 
the paper more self-contained, we recall now the definition of the meaning of the operation 
ffl. To this end, we regard the tessellation Y{t, W) as a frame or primary tessellation and 
associate with each ceU c G Ce\\s{Y{t,W)) an i.i.d. copy Yc{s,W) of Y{s,W). Another 
tessellation Y{t,s, W) of W is now defined as 

Y{t,s,W):^Y{t,W)U [j icnY,{s,W)), 

ceCe\ls{Y{t,W)) 

i.e. we consider the union of Y{t, W) and the cut-outs of Yc{s, Wys within the cells c of 
the primary tessellation (note that in the definition we have used the interpretation of a 
tesseUation as a random closed subset of W). We say that Y{t,s,W) is the iteration of 
Y{t,W) with Y{s,W). The remarkable property of our tessellations constructed by the 
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Figure 1: Realizations of a planar and a spatial stationary and isotropic STIT tessellation 

MNW-process is that the outcome Y{s,t, W) coincides in law with Y{t + s, W), i.e. with 
the continuation of the MNW-construction until time t + s. Thus, it is the same either to 
continue the MNW-construction from t until time t + s or to perform at time t an iteration 
of W) with Y{s,W). 

The local properties established so far, can be extended to the whole space, since the 
random tessellations Y{t, W) are consistent in W, by which we mean that for any fixed 
t > andW CW CR'^ it holds 

Y{t, w')nw = Y{t,W), 

where = stands for equality in distribution. This implies - in view of the consistency 
theorem [11, Thm. 2.3.1] - the existence of the whole-space tessellation Y{t) such that 

Y{t, W) = Y{t) n W for each compact convex W cM. . 

It directly follows from the Markov property that the random tessellations Y{t) are stable 
with respect to the operation ffl, i.e. 

Y{t) =n{Y{t)m...mY{t)), riGN, 

V ' 

n 

where n{-) means the dilation with a factor n. This property explains also the abbrevi- 
ation STIT, because the last equation is a classical probabilistic stability relation. The 
random tessellations Y{t) share another important property, namely that the intersection 
of Y{t) with a j-dimensional plane ii^j, 1 < j < ci — 1, is again an iteration stable random 
tessellation. More precisely, we have 

r(t)nE, =r(7,t,i?,) (i) 
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and the integral-geometric constant 7^ is given by 



., r(ai)r(f) 

r(|)r(^)' ^' 

see Eq. (3.29T) in [7]. Moreover, it is easy to see from the properties of the capacity 
functional of Y{t), compare with [8, Lem. 5(ii)], that STIT tessellations have the following 
scaling property: 

tY{t)^Y{l\ (3) 

i.e. the tessellation Y{t) has, upon rescaling with the factor t, the same distribution as 
y(l), that is the STIT tessellation with surface intensity 1. 

We close this section by mentioning that the random tessellations Y{t) have another im- 
portant property, namely that Y{f}) has Poisson typical cells or Poisson cells for short, see 
Lemma 3 in [8] or part (1) of Proposition 1 below. This is to say, the typical cell of Y{t) 
has the same distribution as the typical cell of a stationary and isotropic Poisson hyper- 
plane tessellation with surface density t (see Chap. 10.3 in [11] for references about the 
classical Poisson hyperplane model). In particular, this fact combined with the intersection 
property from above shows that the typical cell of the lower dimensional STIT tessellation 
Y{t) n Ej, with Ej as above, has the same distribution as the typical cell of a Poisson 
hyperplane tessellation in MP with surface density 'fjt. More formally, we define the cell 
intensity measure M^(*'^) of the STIT tessellation Y{t, W) by 

MnW:=E Yl (4) 

ceCe\\s(Y(t,W}) 

and its (d — l)-dimensional maximal polytope intensity measure W^^^f^^ by 

f2';"'' = e E s„ (5) 

/eMaxPolytopeSrf_i (y (t,VF)) 

where stands for the uni mass Dirac measure concentrated at (•). Moreover, we let 
j^PHT(t,Ty) ^Yie cell intensity measure and F™^*^*'^) be the {d— l)-facc intensity measure 
of a Poisson hyperplane tessellation PHT(t, W) within W G having intensity measure 
tA, which are defined similarly to M^*^*''^^ and W'^^f^^ above (clearly, the (rf— l)-dimensional 
maximal polytopes have to be replaced by the set of {d — l)-faces of PHT(t, W)). These 
definitions bring us in the position to reformulate special cases of Theorems 1 and 2 from 
[12] adapted to our later purposes: 

Proposition 1 We have 

(a) M^(*'^) = mP«^(*'^) and (b) F^i^'^^ = ^F^^^^'^^ds. 

Jo -5 
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2.2 Martingales in the MNW- Construction 

As already noted above, the MNW-construction of iteration stable random tessellations 
Y(t, W) in finite volumes 1^ C R*^ enjoys a Markov property in the continuous time 
parameter t. In our previous work [12] we have used this fact combined with the classical 
theory of martingale problems for pure jump Markov processes to construct a class of 
natural martingales associated to the MNW-process. In this paper we only need a part 
of that theory. To formulate it we let be a bounded and measurable functional on the 
space of {d — l)-polytopes in IR'^ and we denote by MaxPolytopeS(i_i(y) the collection of 
(d — l)-dimensional maximal polytopes of a tessellation Y standing for a generic realization 
of Y{t, W) for some t > 0. Write 

S<^ = S^(F):= Yl (6) 

/6MaxPolytopes^_ ^ ( Y) 

Any hyperplane H G \W] hitting the window W is tessellated by the intersection with Y 
and we denote by Cells(F fl H) the set of all {d — l)-dimensional cells of this tessellation 
and introduce 

MY) f E MKdH). (7) 

It is also convenient to introduce the bar notation for centered versions of these quantities 
withF = r(t,W^),that is to say S<^(F(t, W^)) := Py)) -ES<^(F(t, W)) and likewise 

A^iYit, W)) := A^{Y{t, W)) - EA^{Y{t, W)). With this notation, in view of the results 
developed in [12], we have 

Proposition 2 For bounded and measurable functionals 4> and ip on the space of {d — 1)- 
polytopes in W^, the stochastic processes 

^^{Y{t,W))- [ A^{Y{s,W))ds (8) 
Jo 

and 

f:^iY{t,W))f:4Y{t,W)) - [ A^4Y{s,W))ds 

Jo 

- f [A^{Y{s,W))f:^{Y{s,W)) + A^{Y{s,W))f:^{Y{s,W))]ds (9) 
Jo 

are martingales with respect to the filtration '^t induced by {Y{s, W^))o<s<t- 

Sketch of a Proof: In order to make the paper self-contained, we give here the main 
idea of the proof of Proposition 2, although there is some overlap with [12]. 
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At first, it is a direct consequence of the MNW-construction that the generator L of the 
pure-jump Markov process Y{t, W), i > 0, is given by 

J\W] r^r.'^.r^r.s 



f€Cclls{YnH) 

where Y stands for some instant of Y(t, W) and F is a bounded and measurable function on 
the space of tessellations of W. Applying now the classical Dynkin formula (see for example 
Lemma 19.21 in [4]) with F — gives the martingale property of the random process 

(8). (It should be mentioned that so defined F need not to be bounded in general. However, 
this technical difficulty can be overcome with a suitable localization argument.) 
For the second statement we consider the time-augmented random process {Y[t,W),t), 
t > 0, which has generator L' given by 



L'G{Y,t) = [LG{;t)]{Y) + 



for appropriate functions G{Y, t). Now, using Dynkin's formula now for the product process 
{Y(t,W),t) and with G{Y,t) := (S^(r) - ES^(F))2 = S^F), Y = Y{t,W), gives (again 
after a suitable localization argument) the martingale property of the random process 

T^liYit^W))- f A^.{Y{s,W))ds-2 f A^{Y{s,W))T.^{Y{s,W))ds. 
Jo Jo 

For another functional ijj we can apply the latter property for (p + ijj and (p — ijj and subtract 
the two results, which shows the martingale property of the random process (9). □ 



2.3 Mean Values for Intrinsic Volumes 

In this subsection we discuss certain basic first-order consequences of Proposition 2 to be 
of use for our future reference. We will write Vj for the intrinsic volume of order j with 
j = 0, . . . ,d — 1 (for the standard definition of these functionals we refer to [11] and the 
references cited therein). Further, write 

F,(r) := J2 ^^•(^)- 

ceCells{Y{t,W)) 

Note now that whenever a new facet / splits a cell c into and c" of Y{t, W) giving rise 
to a new tessellation Y', we have 

Fj{Y') - Fj(Y(t, W)) = Vj(c+) + Vj(c-) - Vj(c) = y,-(c+ n C-) = Vj(f), 

since Vj has the valuation property. Consequently, constructing the tessellation Y — 
Y{t, W) by successive cell sphts we easily get 

f,{y) = j:v,{y) + v,{w). (10) 



8 



It is our aim to relate Fj{Y) and Fj+i{Y). For this purpose we use a special case of 
Crofton's formula from classical integral geometry, which reads 

/ V^{KnH)K{dH) = ^^+,V^+,{K), (11) 

J[K] 

where K dW^ is a, convex body and where the constant 7^- is given by (2), see [11, Thm. 
5.1.1]. Applying this formula to Fj{Y) yields 

Ay^{Y) = I V,{f)A{dH)= f Yl V,{cnH)A{dH) 

I'^] /eCeiis(yni?) 1^1 ceCeiis(y) 

= Y 7.+i^.+i(c)=7.+ii^,-+i(l^). (12) 
ceCeiis(y) 

To streamline the notation, we shall write Ey^..^ := Ey^. (F(s, W)) and Avy,s '■= Ay. {Y{s, W)) 
below. Note that, upon combining (10) and (12), it follows from (8) with (j) — Vj that 

^vr,t-lj+i [ ^Vj+v,sds - t^j+iVj+i{W) (13) 
^0 

are all Srniartingales for j = 0, . . . , c? — 1. In particular, bearing in mind that Ey^j^ = 
we conclude that TjVd-i;t — tjiiVdiW) is a S^rinartingale, in particular EEy^_^;t = tVd{W), 
where we have used 7^^ = 1. The latter equation is extended by 

Proposition 3 For j e {0, 1, . . . , 0? — 1} we have 

^^Vr,t-[ n 7.+. 1 Jd^Vd(W)+ Y (117.+^ 1 ^^i+^W 

with the convention that HlLi ■ ■ ■ = 1 and Yl^=i ■ ■ ■ = 0. 
Proof: Taking expectations in (13) we get 

i;t^ f lj+i^^Vj+us,dsi + -fj+itVj+i{W). 
Jo 



EEy. 



Continuing recursively by applying (10) and (12) we end up with 

ft /•Sd-j-2 



EEt 



pt rsd-i-2 

■■■ 7j+l • • • ld-l^^Va-r,sa-j-idSd-j-l ' ' ' dsi 

Jo Jo 



+7,-+i • • • ld-i ^^_._-^y Vd-i{W) + ... + 7,+itV;+i(iy). 
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However, the definitions of Y{s, W) and V^-i imply that EEvd-i;s — sVd{W) for any 
0<s<t. Thus, 



rt rSd-j-2 

■ ■■ 7j+i • • • 7d-i^^Va-usd-j-idsd-j-i ■■■dsi 

Jo Jo 



-i-2 



Ij+i ■ ■ ■ ld^iVd{W) / • • • / Sd-j-idsd-j-1 ■■■dsi 
Jo Jo 

and the result follows immediately. □ 

It is interesting to compare the mean value formula from the last proposition with that 
from [12], Section 3.2. We denote by Tpj(Y(t)) the density of the j-th intrinsic volume of 
the collection of {d — l)-dimensional maximal polytopes of Y{t), i.e. 

for arbitrary compact windows W C with V'd(VF) > 0, see [11, Chap. 4.1]. Using now 
Proposition 3 we get 

because of the homogeneity of intrinsic volumes. On the other hand, we have shown in 
[12] that 

where here and in the sequel we denote hy Hi, < i < d, the volume of the i-dimensional 
unit ball. Indeed, these two values are identical, since 

d\ fnd-A'-Ud _ r(^) / r(f) 1 



r(i) V r(i) 1 



/d-i-j 



n 7.+. 



i=l 



{d-j)\ 
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Reiricirk 1 Our first- order formula from the last proposition contains the extra boundary 

d-j-l / i \ 

correction term | ]^7j+ik 1 TjV^+i(W^) in contrast to our mean value formula (I4) 

i=i \fc=i / 

from [12]. These additional terms come from the fact that we consider maximal polytopes 
possibly chopped off by the boundary of the domain W rather than the points of an asso- 
ciated center function for full facets in M.'^. Thus, it may happen for instance that in two 
neighbouring regions one observes two facets which can coalesce when putting these regions 
together into one volume, whence the lower-order finite volume corrections arise. 



3 Second Order Characteristics for Intrinsic Volumes 

In this Section we develop a full second-order theory for intrinsic volumes of stationary 
and isotropic random STIT tessellations. This is first done exactly in Subsection 3.1 for 
Y{t, W), t > inside a bounded convex observation window W C M.'^. Then, in Subsection 
3.2 we derive the corresponding asymptotic expressions for Y{t, Wr), Wr :— RW, i? — > 00. 



3.1 Exact Expressions 

To proceed with second-order calculations we shall use the notation already introduced in 
Subsection 2.3 above. Observe first that in view of (10) the relation (12) simphfies and we 
obtain 

A,.(F) = 7,+iSv,^,(F). (15) 

Recalling that Sy^.^ := Hy^(Y{s,W)) and Ay^.^ := Avj{Y{s. W)), putting likewise Ay^Vjis '■— 
AviVj{y{s, W)) and using (9) in Proposition 2 with 4> = Vi and ip = Vj, we see that 

f^Vi;t^Vj;t - / AviV^^sds- / [Avi;s^Vr,s + Avj;s^Vi;s]ds 

Jo Jo 

is a martingale with respect to the filtration 9=f induced by {Y{t,W))Q<s<t- Substituting 
(15), this implies that 

_ _ /■* /■* 

Sy..+Sy . 



/ AviVj;sds- / [^i+iT,Vi+j^-s'^Vj;s + lj+l'^Vi;s'^Vj+v,s]ds 

Jo Jo 



is an Q^^-martingale as well. This is a crucial formula because, upon taking expectations, 
it allows us to express Cov(Sy-.t, Sy .() in terms of corresponding covariances with indices 
{i,j -\- 1) and (j, i -|- 1). In other words, we get 

Cov(Ey..j,Ey..t) = / EAviVr,sds+ / Cov(Sy.^^.„ Sy..,) + 7^+1 Cov(Sy,.„ Ey.^j.,)]cis. 

^0 Jo (Ig) 



11 



It is important to observe that this recursion terminates because = 0, which allows 
us to provide an explicit expression for the covariances Cov(Ey..t, Ey^.-t). To obtain the 
desired formula, denote 

rt rsi rsn-i i rt 

T^if; t):= / . . . / ...ds, = — / {t - sY-'f{s)ds 

Jo Jo Jo \1^~^)-Jo (^lY'j 

for any / : [0, t] — )■ M and n eN for which the iterated integral exists. Then we claim that 
for /c, / e {0, . . . ,d-l} it holds 



(^r) (n^'^-^) (ri^'^-)^'^'(Var(E^,_,(.));^)+ 



0<m<fc. 0<ri<i ^ ^ \i=m+l / \j=n+l / (J-Oj 

(m,n)^(0,0) 

with the convention that (g) = 1 and Ylt+i ■ ■ ■ — Y[\+i • • • = ^- Whereas this can be 
readily verified by a straightforward induction in view of (16), there is a more natural 
way to see it. In fact, in the course of recursive applications of (16) the covariance Ck^i '■— 
Cov(Sy^_^_^, T,Vd-i-i) is represented in terms of Ck-i,i and Ck,i-i which can be interpreted in 
terms of lattice walks on pairs of indices from (A;, I) to (0, 0) where only steps (i, j) — >■ (i — 
and (i, j) — > — 1) are allowed, each step {i,j) — > (i — involving multiplication 
by 7rf_j plus integral iteration and each step — )■ — resulting in multiplication by 
7rf_j plus integral iteration. There are C^^') such walks, whence the first line in (18) follows. 
The second line of (18) is due to additional [!^E.Av, , , ,v, . , ,-sCis terms in (16), which 
are born at all times (m', n') of the lattice walk discussed above, with m' standing for the 
number of — >■ (i — 1, j) steps and n' for the number of — >■ — 1) steps. Note 
that no additional term is born at (m' = k, n' = I) though, as it corresponds to having 
reached the (co)variance Cov{Hy^_^, Ey^ J = Var(Ey^_J, which does not get expanded any 
further. In (18) we have substituted m = k — m' and n = l — n' for notational convenience. 
Applying now (16) ioi i = j = d — 1 and using that Ey^.j = we see that 

Var(Ey,_,;,)=X^(EAy^._^^(.);i), (19) 

whence we finally get from (18) the following exact variance expression: 

Theorem 1 The covariance between the intrinsic volume processes Ey^_j_^.i and '^Vd_i_i;t 
for k, I E {0, . . . , d — 1} of a stationary and isotropic random STIT tessellation Y{t, W) is 
given by 
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tti'^liZ'") i n ^-.) f ri 7.-.)i'+'-"-»«(E^v._._„.._._„(.,;*) 

m=0 n=0 ^ ' \j=m+l / \j=n+l / 

with the usual convention that (q) — 1 anc? Hifc+i ■ ■ ■ — Yll+i ■ ■ ■ = 1- 

Unfortunately, wc arc not able to make the covariance in the above theorem any more 
explicit. This is due to the presence of mixed moments E,Ava_-^_^Vd-i^„;t, rn,n> 0, whose 
evaluation is technically related to the problem of providing a general formula for joint 
moments of lower-order intrinsic volumes oi {d — l)-dimensional Poisson cells (arising as 
hyper-planar sections oiY(t, W), which are in addition possibly chopped off by the bound- 
ary of the window W) which is not currently available in required generality up to the best 
of our knowledge. Fortunately though, the offending mixed moments turn out to be of 
negligible order in large window size asymptotics and we are able to provide fully explicit 
asymptotic formulae in Subsection 3.2 below. 

We would like to point out that in the special case d — 2 our Theorem 1 reduces to The- 
orem 1 and Corollary 1 in [13] where we have studied variances and central limit theory 
for maximal edge (I-segment) lengths and vertex counts and where we could give fully 
explicit exact variance formulae exploiting the particular features of the planar setting. In 
Corollary 2 ibidem we also have provided asymptotic expressions for these variances for 
sequences of growing observation windows. For this reason, in our asymptotic considera- 
tions below we will restrict to the case d > 2 where the general asymptotic expressions are 
essentially different from those arising in the exceptional planar case. This is due to the 
variance dichotomy established in [12] and mentioned in Section 1 above. 



3.2 Asymptotic Expressions for d > 2 

Let W G M.'^ he a compact convex window and consider the sequence Wr := RW. We 
will write from now on ^Jf instead of A^{Y{t, Wr)) and likewise Sjf for E^(F(t, Wr)) 
in order to emphasize the dependence on R. Our main interest in this section is to derive 
from Theorem 1 an asymptotic expression for the covariances Cov(S^^^_^.^, S|^^^_^.J as 
it! — >■ oo. We start with the following 

Proposition 4 For k,l,m,n & {0, . . . ,d — 1}, m < k, n < I, with t fixed we have 

r 0(i?2(««-i)-™-«), ifm + n<d-3, 
Xk-,l-m-n^l^^^Wn^ ^v^, , „.,^.yt) = I 0{RHogR), m + n = d - 2, 

( 0{R'^), ifm + n>d-l. 

Proof: Observe first that, in view of (17), the integral x'^+'~™~"+^ only involves integra- 
tion with respect to the variable s and does not affect the order in R. Thus, it is sufficient 
to prove the Proposition for A; = m and / = n which we shall henceforth assume without 
loss of generality. We have 
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where 



4'l= f y .JS 

Jl R 



'1/R 

and 

/ EA^'' y .As. 
Jo 

To provide a bound for we need some additional notation. Write w,n;s for the common 
value of 

Eyd_i_„(TypicalCell(y(s) n //))Ki_i_„(TypicalCell(y(s) n H)) 

with H ranging through hypcrplanes in R*^ and where TypicalCell(y (s) fl H) the typical 
cell of the sectional STIT tessellation Y{s) fl H. Using the scahng property (3) and the 
homogeneity of the intrinsic volumes we readily get 

^ ^ m+n-2(d-l) /2n^ 

Sm,n;s Sm,n;l- \ ^/ 

Indeed, 

<im,n;s = El^rf_i_^ (TypicalCell(r (s) n i/))Vd_i_„(TypicalCell(r (s) n H)) 

= EVd-i-m TypicalCell(sr (s) n H)^ Ki-i-n iypicalCell(sr (s) n H) 

d—l—m+d—l—n 

EFd_i_^(TypicalCell(r(l) n i/))l^rf_i_„(TypicalCell(F(l) n H)) 

— * Sm,n;l- 

To proceed, write, recalling (7), 

Wb\ /6Cells(y (s,Wfl)ni?) 



/GCeiis(y(s)n/r), /nw«^0 



= \ ^^^^E Vold-i(TypicalCell(F(s) n i/) © {-{Wr n H)))A{dH), 

J[Wr] ^0,d-l;s 

where © stands for the usual Minkowski addition. Recalling that TypicalCell(F(s) n H) 
is Poisson, see Section 2.1, we readily get for s > l/R, 

E Vold_i(TypicalCeU(r(s) n H) ® {-{Wr n H))) = 0(i?'^-^) 

and thus we conclude in view of (20) that, for s > 1/i?, 
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Consequently, 



= o{r^ ) = I OiR'^log R), if m + n = d - 2, 

^ "^VK / ( 0(i?''), ifm + n>d-l. (21) 

(2) 

Next, we find a bound for 3)^.^. To this end, we note that during the time interval [0, l/R] 
there are 0(1) cell splits within Wr in the course of the MNW-construction and hence the 
expectations of sums 

l^d-l-m(/)K«-l-n(/) 

feCens{Y{s,WR)nH) 

are of order 0{Vd-i-m{WR)Vd-i-n{WR)) = 0{R^^'^-^^-") uniformly in hyperplanes H e 
\Wr] and in s e [0, l/R]. Thus, using that A([W^i?]) = 0{R) and recalling the definition (7) 
we obtain 

5g ^oi^- R^(d-l)-m-n . ^ ^ ^ ^^^2(d-l)-m-ny ^22) 

Putting (21) and (22) together completes the proof of Proposition. □ 

This implies that asymptotically, as it! — >■ oo, the terms appearing in Theorem 1 with 
n, m > are of order at most 0{R^'^~^) and thus neghgible compared with the leading 
7^'^(fi-i)_tcrm. For this reason, the asymptotic behaviour of Cov(Ev'^_^_^;t, Ey^ -^ j^^) is domi- 
nated by the term with m = n = 0, which is of order 0(/?^*^^~^)), where by ©(■) = 0(-)nr2(-) 
we mean quantities bounded both from below and above by multiplicities of the argument 
in the usual Landau notation. For this dominating case Proposition 4 is refined by 

Proposition 5 We have for /c, Z e {0, . . . , c? — 1} 

1 (A; + 0! 2 ^ ^ 

Proof: First, recall from Thm. 4 in [12] that the variance of the total surface area of 
Y{t, W) equals 

d — 1 f f 1 — e 
Var(Ey,_,;,) ^ dxdy. (23) 

^ JwJw \\x-y\\ 



I [ — ^dxdy + OiR^"^-^). 

Jw Jw \\x - y\\ 
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The main argument for (23) reads as follows: At first, we use the fact that STIT tessella- 
tions have Poisson typical cells and find 

Ay.jY{s,W)) = I Vl,{f)A{dH) 

-^[W] f(zCens(Y(s,W)nH) 

y are in the same cell of Y{s, W) fl H]dxdyA{dH) 



[W] JwnH JwnH 

e-'H[m)dxdyA{dH), 



[W] JwnH JwnH 

where ![•] stands for the usual indicator function and xy for the line segment connecting 
X and y. Using now (9) with (j) — ip, taking expectations and noting that the mixed term 
vanishes leads to 

Va.r{Ev,_,;t)= f Ay2 (Y{s,W))ds = fl I f e-'''^^'^^^dxdyA{dH) 
Jo Jo J \w] JwnH JwnH 



'0 J[W] JwnH JwnH 

f f r 1 — e~*^([^l) 

= / / / — . — dxdyA{dH). 

J [W] JwnH JwnH A ( [xy\ ) 

The latter integral can be transformed into (23) by using an integral-geometric formula of 
Blaschke-Petkantschin-type and the Crofton formula (11) - see [12] for details. 
Recall now (19) and write 

l'^-^^+\EA^f .^.);t) =X^+'(Var(S;^«^^(.P;t) = i?^('^-^)X'^+'(Var(SVl,^(.)«); t), 

(24) 

where the last equality follows by the scaling property (3) of STIT tessellations. Thus, 



i?2(.-i) d-1 f f 



-Rs\\x-y\\ 



^'^'^'i^Kf i-yt) = / / {t-s)'^^-' dsdxdy. 

''d-i'U + 2 JwJwJo \\x-y\\ 

Now, the binomial theorem implies that there exists some constant c > such that 



Jo 



^ \ — Q-c\\x-y\\Rs 



\x - y\\ 



ds 



1 t^+^c^+i \\x - E}'+^ + ce-'^b-fll*^ + 0(i?*^+'-^) 
{k + l) \\x-y\\''^^^\^+^R^+i 



{k + l) \\x-y 



|2 



+ 0{l/R) 



with c = In view of (24) this proves the desired result. □ 

Consequently, by combining Propositions 4 and 5 with Theorem 1 we get 
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Corollary 1 Asymptotically, as R ^ oo we have for k,l & {0, . . . ,d — 1} 

1 d-1 



C*')(nH(ri^(*+o! 2 

where E2{W) denotes the 2-energy ofW, i.e. 



-t''+^R'^^-^^E2{W) + 0{R 



2d-3\ 



E^iW)^ [ [ —^—^dxdy. 
Jw Jw \\x — y\\ 



y\ 

The affine Blaschke-Petkantschin Formula [11, Thm. 7.2.7] can be used to provide an 
integral-geometric expression for the measure-geometric energy functional E2. In fact, we 
have 

where Id-i(W) denotes the {d — l)-st chord power integral of the convex body W in the 
sense of [11], page 363, i.e. 



where C denotes the space of lines in R.^ with invariant measure dL and is the volume 

of the (i-dimensional unit ball. This means that in the asymptotic covariance formula 
from Corollary 1 the dependence on the geometry of W is encoded by the non-additive 
E2{W) or equivalently by Id,-i{W). We will from now on use the representation in terms 
of chord power integrals, as it allows an easier comparison with other tessellation models. 
Summarizing, this yields the following 

Corollary 2 Asymptotically, as R ^ 00, we have for k,l & {0, . . . ,d — 1} 



k \ / I 



c-(e;£._,,. = ^ (117.-. ) ( n^-. ) ^''-^m^"'"-'' + 

and for k — I 



\li Ul / kill 

^1=1 / \j=i 



2 

^2k 



Especially for the practically relevant case d = 3, we have for A; = 1 and k — 2 



Var(S;;-) = -t%{W)R' + 0{R% 
V^^(^S) = ^t%{W)R' + 0{R'). 
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In general, Id-i{W) or equivalently E2{W) is rather difficult to evaluate explicitly. How- 
ever, for the unit ball W — B'^ in M!^ we have by applying [11], Theorem 8.6.6 (with a 
corrected constant) 

For example for d = 3 this gives us the value hiB^) — 47r^. For the interpretation 
of computer simulations it is of particular interest to evaluate the chord power integral 
liiCa) for a 3-dimensional cube with side length a > 0. It can be shown that the 
numerical value of hiC^) is given by 5.6337 • a'^. 



4 Second Order Structure of Face Measures 

In this section we focus our attention on the spatial pair-correlation structure for the 
processes of maximal polytopes of arbitrary dimensionalities induced by Y{t), arising for 
J = 0, . . . , c? — 1 as random j- volume measures concentrated on the union of all j-faces of 
{d— l)-dimensional maximal polytopes of Y{t). The nature of the so defined face measures 
is somewhat different than that of cumulative intrinsic volume processes considered above, 
in particular the face measures keep track not just of the cumulative numeric characteristics 
of STIT tessellations but also of their spatial profile, moreover it should be emphasized 
that in general the total mass of a face measure may be quite unrelated to the cumulative 
intrinsic volume of the corresponding order, as for example Svo(y(t, W^)), the number of 
{d— l)-dimensional maximal polytopes of Y{t) inW G M"^, is not deterministically related 
to the number of vertices of Y{t, W) as soon as d > 2 etc. However, we decided to consider 
the face measures in this paper as they are of interest in their own right and supplement our 
results from the other sections, showing the power and versatility of the methods developed 
in this paper. 

For J = 0, 1, . . . , d — 1 we consider the (random) j-th order face measure 
V,,:= Yl Vol,(-ne)= Yl ^/ 

eeMaxPolytopes^ {Y{t)) /eMaxPolytopes^j.i (Y (t)) 

where Volj(-ne) is the e-truncated j-volume measure [Volj(-ne)](74) = Volj(Ane), A C R'^, 
whereas 

vf:= Y Vol,(-ne), 

egFaccsj (/) 

where by FaceSj(/) we mean the collection of all j-faces of the (rf— l)-dimensional polytope 
/. We also abuse the notation by putting Vd-t '■= Vold(-)- We shall be interested in the 
covariance measures Cov(Vi;t, Vj-t), ^, J = 0, . . . , d — 1, on (M*^)^ given by 

{g ® h, Cov(V,;i, V,.t)) = CoY{{g, Vr,t), {h, Vj-t)) 

for all bounded measurable g,h : ^ with bounded support, where the standard 
duality notation (•, •) is used for integration (0, jj) — j 4>diji. Denote by Q^i^ the law of the 
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typical {d — 1) -dimensional maximal polytope, of the STIT tessellation Y{t) and let ^Y{t)^ 
be the corresponding facet density. The following Proposition is crucial for this section: 

Proposition 6 For i, j e {0, . . . ,d — 1} we have 

[ [7i+i Cov(H+i;„ V,-;,) + 7,-+i Cov(H;„ V,+i;,)]ds (25) 
with 7j+i and 7^+1 given as in (2). 

Proof: For bounded measurable and boundedly supported g and and for i, j G {0, . . . , d— 
1} consider the facet functionals 

J!{f)^{9M)- E (^,Vol,(-ne))= E I g{x)Vol,{dx) 

eeFacesi(/) eeFacesi(/) ^ 

and defined analogously. Then, choosing some compact convex W containing the 

supports of g and h in its interior we see that, recalling (6) and (7), 

Aj.{Y) = I E JUmidH) = E L [ Voki- n (e n H))A{dH)\ 

>^[M^] /gCells(yn//) eeMaxPolytopeSi+i(y) \ l"^] ' 

= 7,+i E (5,Vol,+i(-ne)) = 7,+iE^^^(r), (26) 

eGMaxPolytopeSj.,.! (Y) 

where we have used the Crofton formula (11). Applying now (9) with = J? and ip — Jj, 
we get upon taking expectations 

Coy{Ej.{Y{t,W)),Eju{Y{t,W))) = [ EAj.ju{Y{s,W))ds+ 

[Coy{Aj.{Y{s, W)), Ej,{Y{s, W))) + Coy{Ej.{Y{s, W)),Aj,{Y{s, WMds 
and thus, in view of (26), 

GoY{llj.{Y{t,W)),Y.jH{Y{t,W))) = f\Aj.jH{Y{s,W))ds+ 

^ Jo ^ ^ 

[7m Cov(S^. (F(5, W)), S^.(F(s, 1^))) +7i+i Cov(S^.(F(5, W)), ^j.jY{s, W)))]ds. 
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(28) 



Using that 'Ej?{Y{t, W)) — {g, Vi-t) and similar relationships, we end up with 

{g®h,Cov{Vi.,t,Vj.,t)) = [\Aj.j,(Y(s,W))ds+ 

Jo ' ^ 

(9®^^ bi+i Cov(Vi+i;„ Vj-s) + 7i+i Cov(Vi;„ Vj+v,s)]ds^ . (27) 
Thus, to estabhsh (25) it is enough to show that 

j\Aj.ju{Y{s,W))ds^(g®K\^i-^ j^J ^P^)®^;y^^)Q2?(d/)dx 

To prove (28) use (1) and write, recalling (7), 

¥.Aj.jn{Y{s,W)) = /" E JfJ^imdH) 

^ feCc\\s{Y{s,W)f\H) 

-II J/Jt(/)M^(^'^)"^(ci/)A(cii/), (29) 
J\w] J 

where M^*^*''^)^'^ is the cell intensity measure for the sectional STIT tessellation F(s, PF)n 
H, compare with (4). Using Proposition 1 (a) we are led to 

I Jf jj^(/)M^(^'^)"^(d/)A(di/) = / / Jf jj^(/)MP«^(^'^)"^(d/)A(d//). 

[W] J J[W] J 

Now, using Slivnyak's Theorem (cf. [11, Thm. 3.2.5]) we obtain for any bounded measur- 
able function on the space of {d — l)-dimensional polytopes 



0(/)MPHT(i,H^)n//(^y)A(rfi7) = / / 0(cni/)A(di/)MP«T(i'^)(dc) 

J J[c\ 



where M^^^^^'^)"^ is the cell intensity measure of PHT(1, W) n H, whereas M^^^^^'^) is 
the cell intensity measure and F^^^*^-*^'*^) is the {d — l)-face intensity measure of a Poisson 
hyperplane tessellation PHT(1, W) within W having intensity measure A. Thus, replacing 
A by sA, < s < t, we arrive at 

/ / <t>{f)yf^^^^'^^-"{dmdH) = 1 / 0(/)Ff 

J[W] J ^ J 

whence continuing (29) we find 

EAj.j.{Y{s, W)) = -J J!J'^{mT,^'''^\df). (30) 



20 



Hence, by applying Proposition 1 (b) we get from (30), 

j\Aj.j.{Y{s,W))ds^ J J!J^{ml^T\df\ (31) 

where, recall, fJI*^'^'* is the [d — l)-dimensional maximal polytope intensity measure for 
Y{t^W), see (5). In view of (31) the required relation (28) follows now directly by the 
definition of typical {d — l)-dimensional maximal polytope [11, Chap. 4.1, Chap. 10] and 
the definition of Jf , upon taking into account that the supports of both g and h are con- 
tained in the interior of W. Putting (27) and (28) together yields (25) and thus completes 
the proof. □ 

To proceed note that Vd-^s is a constant measure and therefore any covariances involving it 
vanish. Thus, using (25) in Proposition 6 and arguing as in the derivation of Theorem 1 
from the crucial relation (16) we obtain 

Theorem 2 For /c, / G {0, . . . , d — 1} we have 

cov(v._._«,v._._„)= ttf*^!:^"") ( n ^.'-) f ri -^A 

m=0 n=0 ^ ^ \i=m+l / \j=n+l / 



X £ / bSl ® et-ll^'^^'-'"-" {^^i-I^Ql^ht) {df)dx (32) 

with the usual convention that (g) = 1 and Ylt+i ■ ■ ■ — Y[\+i ■ ■ • = 1- 
It is interesting to note that for the particular case A; = Z = we get 
Corollary 3 We have 

Cov(V._,,) = £^ J [vlr^ ® vlrWa^'hdf)dx. 

This means that foik — l — O the covariance measure of the surface area process Vd-i-t 
coincides with that of the surface area process induced by the homogeneous and isotropic 
Boolean model with grain distribution Qji*i and with grain density )^Y{t) '' ^ ^ result first 
estabhshed by Weiss, Ohser and Nagel [14] in the special planar case by completely different 
methods. Recall further from Theorem 3 in [12] that 

where Q^-i stands for the distribution of the typical facet of the Poisson hyperplane 
tessellation PHT(s) with surface intensity s (this is to say, the hyperplane process has 
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intensity measure sA) . But this random polytope has the same distribution Q^JcaiCeii ^ 
as the typical cell of a Poisson hyperplane tessellation in R'^"-'^ with surface intensity ^a-is 
with 7(i_i given by (2). Thus, 

and the covariance measure is reduced to known quantities. 

For k + I > the situation becomes more complicated and mixtures of typical cell distri- 
butions corresponding to different time moments arise in the right hand side of (32). As in 
Section 3 above, also here explicit calculations are precluded for /c + / > due to the lack 
of known formulae for mixed moments of general order intrinsic volumes of Poisson cells. 



5 Central Limit Theory 

In this section we present a central limit theory for the suitably rescaled intrinsic volume 

processes Sy^.^. In strong contrast to the two-dimensional case considered in [13], where a 
classical Gaussian limit behaviour is observed, for the case d > 3 in focus of this paper 
the situation is very different and non-Gaussian limits arise. In this context, to proceed 
with a full discussion below, we recapitulate first some facts already known from [12] in 
a way specialized for our present purposes. Define the rescaled intrinsic volume processes 
(»5^';r)te[o,i] for i = 0, . . . , d - 1 by 

S^^ = R-^'-'^I^Zw/R = R-^'-'^^vAYit + logi?/i?, Wn)l t e [0, 1]. 

Note that the shift by log R/R in time argument here is of technical importance as placing 
the time origin just after the very initial big bang phase of the MNW-construction, where 
the dominating fluctuations of f^Vd-i;{-) arise, but where nothing of asymptotic significance 
happens for ^Vi;{-), i < d — 1. Since the big bang phase evolution for Sy^_^.(.) has been 
considered in full detail in [12, Sec. 5], in this paper we only focus on the later phase 
[log it!/i?, 1]. Putting together the present Section with [12, Sec. 5] reveals remarkable 
richness of the complete asymptotic picture for STIT tessellations in large windows. 
We begin with the observation that follows by the theory developed in Subsection 5.3 of 
[12]: 

Proposition 7 The process (iSyj ^.^)te[o,i] converges in law, as R ^ oo, in the space ©[0, 1] 
of right- continuous functions with left-hand limits (cddldg) on [0, 1] endowed with the usual 
Skorokhod Ji-topology [3, Chap. VI. 1] to the constant process t ^ where ^ :— E{W) 
is a certain non-Gaussian square-integrable random variable with variance ^^E2{W) — 

Proof: The relation (85) in [12] implies that converges in law to ^ as 

R oo. On the other hand, recall that Sy^_^.^ is a martingale in view of (13). Moreover, 



22 



by (23) and the definition of <5^'^.t combined with the scahng property of STIT tessellations 
we have 

Var(5g^,) = YB.r{R-^'-'my^^^{Y{t + \ogR/R,Wn))) 
= Var(Sy,_,(r(i?(t + logi?/i?),l^))) 

-[Rt+\ogR]\\x-y\\ 



2'=d-l I 



1 [ [ 1 — e "^"d , , 

dxdy 



and hence Var(5yj^.j^) — Var(5yj^.Q) tends to as — )• oo. Consequently, the asymp- 
totic constancy of the limit process follows now by Doob's L^-maximal inequality [5, Thm 
3.8(iv)] which completes our argument. □ 



Remark 2 The convergence in law and non-Gaussianity results referred to above come 
from the paper [12]. In order to keep the present work formally self-contained we give here 
a very brief sketch of the main arguments. 

The convergence is readily guaranteed by the martingale convergence theory as ibidem. The 
crucial point is the non-Gaussianity. Here, the idea, as described in detail in Section 5.3 
of [12], relies on showing that ^ has its tails much heavier than normal random variables. 
We proceed, roughly speaking, by constructing a suitable class of initial frame tessellations 
of the body W with the properties that 

• the number of frame hyperplanes hitting W equals N eN, 

• the frame hyperplanes intersect within W only very seldom, which implies that W 
gets subdivided into Q{N) cells (recall the Landau notation). 

This can be achieved by choosing d principal directions and keeping all hyperplanes in the 
collection under construction approximately aligned to these directions. Importantly, it 
turns out that this way we can get the lower bound ex.p{—0{N log N)) for the probability 
that the real initial frame arising in the very initial phase for Y{-, W) - usually called big 
bang, see below - upon suitable spatial re-scaling does fall into the desired class. Now, 
denoting the last event by S^, we see that on £m the total deviation ^ a; Sy^^.^ gets 
decomposed into 

• a sum of independent contributions coming from all respective cells Wi of the frame 
tessellation, 

• a correction term of order Q{N) due to the centred joint contribution of all hy- 
perplanes of the frame tessellation (observe that each such hyperplane necessarily 
contributes 0(1), whereas = 0(1)/ 
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Finally we can use the knowledge o/Var(^) anc? Var(^j) and properly exploit the indepen- 
dence of the random variables to conclude that on S^- we can find a further sub- event 
S'j^ C gj^ with F{S'^) = e(P(^7v)) = exp(-0(A^logA^)) such that > N} D 8'^ and 
hence > N) = exp(— 0(A^log iV)), which could not hold if ^ were Gaussian. We refer 
the reader to [12] for further details. 

Having characterized the asymptotic behaviour of {Sy'^^.-i)te[o,i\ we are now prepared to 
describe the full joint asymptotics of all intrinsic volume processes, which is the main result 
of this section. 



Theorem 3 The vector [Sy'J^^.^, Sy'J^^_^.^, . . . ,Sy^)te[o,i\ converges in law, as R ^ oo, 
the space "^([0, 1]; W^) of -valued cddldg functions on [0, 1] endowed with the usual Sko- 
rokhod Ji-topology [3, Chap. VI. 1], to the stochastic process 



m 



t ^ ( C,7<i-i*C, 2! ^' 3! ^' ■ ■ ■ ' — {d-l)\ — ^ ^ ' ^ ' ^' 



(33) 



Before proceeding with the proof we discuss certain striking features of the phenomenon 
described in Theorem 3. Namely, although all intrinsic volume processes 5^.'^ exhibit 
fluctuations of order R^"^, the nature of these fluctuations differs very much between 
i — d — 1 and i < d — 1. 

• As shown in Subsection 5.3 in [12] and in Proposition 7 above, the leading-order 

deviations of Sy^^.^ in large R asymptotics arise very early in the course of the 
MNW-construction, in its initial stages usually referred to as the big bang phase. 
Here, this is the time period [0, logi?/i?]. During the later stages of the construction, 
i.e. say the time interval (log R/ R, 1], the variance increase is of lower order and any 
newly arising fluctuations are negligible compared to those originating from the big 
bang. In the asymptotic picture this means that the initial fluctuation remains frozen 
throughout the rest of the dynamics, whence the constant limit for (»5^'^.^)tg[o,i] as 
i? ^ oo. 

• In contrast, the leading-order deviations of Sy.'.^ for i < d — 1 arise and cumulate 
constantly in t with deterministic polynomial rates depending on i, times the initial 
big bang fluctuation of the process E^^'^.^. -j which, in this sense, stores the entire ran- 
domness of the intrinsic volume vector. The mechanism determining the dependence 
of fluctuations of intrinsic volume processes of orders i < d — 1 given those for d — 1 
and the resulting form of the limit process (33) will be discussed in the sequel. 

Proof of Theorem 3: The crucial step of the proof relies on considering for each j = 
0, . . . ,d — 1 the auxiliary process 

T,vj;t ■= ^Vj;t - / Avj-sds, (34) 
Jo 
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which is a centred $5t-martingale by (8) and which is the same as 

tvj;t = - 7^+1 / f^v^^usds (35) 



by (15). The idea is to show that for j < d — 1 the processes f^Vj;t are of neghgible order 
in large R asymptotics. Indeed, upon squaring and taking expectations we get 



E(Ey^;t)^ = Var(Ey^;t) - 27j+i / ET,Vj-t^Vj+i;sds + 'y]+i / EEvj+^-s^Vj+uududs. 

Jo Jo Jo (36) 

Using that 

E(Sy.;t|Sjs) = ^Vf,s + ^{^Vi+y,u\'^s)du, 

as follows by the martingale property of f^Vj;ti and inserting this to (36) we are led to 
E(Ev.,i)' = Var(Ey..i)-27,+i / EEy..,Ev..^,.,ds- 

^0 



27J+1 / j E(E(Ey.^^.„|$>,))Ey.^,.,duds + 7j+i / / ^'^Vi+v,s^Vi+v,ududs 

Jo Js Jo Jo 

with the last two terms cancelling. Thus, 

E(Ey.;t)^ = Var(Ey.;t) - 27^+1 / Cov(Ev/.;3, Ey._^^;s)rfs, 

Jo 

whence, by (16), 

E(Ey..t)2= !\Ay..As. (37) 







In view of Proposition 4 for j < 0? — 1 we have 



and thus the relation (37) implies 

Var(E^«)=o(i?2('^-^)), te[0,l] 

for J < d — 1. Hence finally 



Esup(S|:r«)2 = o(i?2('^-^)) (38) 

tG[0,l] 
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for j < d — 1 by Doob's L^-maximal inequality [5, Thm 3.8(iv)] applied to the martingale 
defined as in (34, 35) with W replaced by Wr there according to our usual notational 
convention. 

With (38) it is now easy to complete the proof. Indeed, since the normalization in the 
definition of intrinsic volume processes Sy^^ involves a prefactor R~^'^'~^\ the relation (38) 
allows us to recursively substitute 

7i+i / ^?-Z;s^^ 
Jo 

for Sy^^ , as soon as j < d—1, without affecting the large i? — )■ oo asymptotics in law (note 
that the technically motivated shift by log R/R — ?> in time argument of the rescaled pro- 
cess Sy'^ is asymptotically negligible as inducing only a negligible L^-difference precisely 
calculated in Theorem 1). The application of such recursive substitutions combined with 
Proposition 7 completes the proof of Theorem 3. □ 

It is crucial to emphasize at this point that, in intuitive terms, under the normalization 
of Sy'^ , j = 0, . . . ,d — 1, the mechanism governing the rise of fluctuations of intrinsic 
volume processes of orders j < d — 1 given those for V^-i reduces effectively to the simple 
approximation 

J 

and its recursive application. Of course this simple approximation follows itself by rather 
non-trivial arguments above. Clearly, this description only characterizes the leading order 
fluctuations as considered in Theorem 3 although our tools should allow a more delicate 
characterization of non-leading lower order fluctuations as well, see e.g. Subsection 5.1 in 
[12] for the particular case of Vd-i. 
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